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We study the time development of correlation functions at both zero and 

finite temperature with Shibata-Hashitsume's projection operator method 
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, and carry out the renormalization of ultraviolet divergence that appears in 



a time-dependent frequency shift, using a mass counter term. A harmless 
divergence of log i-type remains at an initial time t = after the lowest order 
renormalization. 

11.10.Gh,03.70.k,11.10.-z,05.70.Ln 



Calculating various physical quantities in a non-equilibrium system is a very interesting 
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■ problem. However, it is well-known that there are various problems in calculating time 
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dependent physical quantities, e.g., "the surface divergences" |I[ and the short time behavior 
of an unstable system which is related to "Quantum Zeno's paradox" j2|. In these problems, 
one encounters peculiar divergences which do not appear in usual S-matrix theory where 
the relevant time interval is [—00,00]. To avoid these divergences, for example, diffuse time 



^ ■ boundaries have been introduced |] @. 

In this paper, we show first that a divergence appears at the initial time t = in 
correlation functions obtained by using Shibata-Hashitsume's projection operator method 
0- and then show that it becomes harmless in calculating some physical quantities in 
a super-renormalizable model by mass renormalization. Shibata-Hashitsume's projection 
operator method is one of various formulations of non-equilibrium statistical mechanics: 
e.g. linear response theory M, Caldeira-Leggett's theory [[7), closed time path method M 
0, projection operator method JlO[] — ||, and so on |13f jl4l . 



Advantages of Shibata-Hashitsume's projection operator method in our analysis are that 
its generic formulation is independent of the explicit form of the projection operator and 
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a non-perturbative effect can be taken into account by solving the Heisenberg equation of 
motion obtained by a projection. 

Our strategy is to study quantum field theory as an initial value problem, that is, we 
prepare an initial state and determine the time evolution by solving the Heisenberg equation 

jO(t)=i[H,0(t)] (!) 
= iLO(t), (2) 

where L is the Liouville operator. The Heisenberg equation has complete information on 
the time development of the operator, but is difficult to solve exactly. So it is necessary to 
make some approximations. For this purpose, we introduce generic projection operators P 
and Q which have the following general properties: 

P 2 = P, (3) 
Q = l-P, (4) 
PQ = QP = 0. (5) 

By these projection operators, we can derive the following two equations. 

d 

di 
d 



—e iLt P = e iLt PiLP + e lLt QiLP, (6) 
at 



dt e iLt Q = e iLt PiLQ + e iLt QiLQ. (7) 

We can solve Eq.(^) for e tLt Q and substitute the solution into Eq.(^). Then we can transform 
the Heisenberg equation ([I]) into the following form: 

jO(t) = (e iLt P + g e *«*)_L_iLO(o), (8) 

where 

E(t)= f ' dre- iL{t ~ T) PiLQe lLQ{t ~ T) . (9) 
Jo 

Now the total Hamiltonian can be divided into three parts: system(s) which we are 
interested in, environment (e) and interaction(se) between system and environment: 



H = H S + H e + H se . (10) 

Let us specify the projection operator. We are interested in the detailed behavior of the 
system, but not in that of the environment. Therefore, we assume that 

P = PsPe, (11) 

where p,p s and p e are the density matrices of the total system, the system which we are 
interested in and the environment, all given at t — 0, respectively, and then define the 
projection operator as 

PO = Tr[ Pe O] = (0) e (12) 

for any operator O. By this projection, we can replace an operator of the environment with a 
c-number. Furthermore, we consider, for simplicity, a special case where p e has only diagonal 
components when it is expanded in terms of the eigenstates of H e . Using these properties 
of the projection operator, we have the following relations which play an important role in 
the following calculations: 

PL S = L S P, (13) 
PL e = L e P = 0, (14) 

where L a O = [H a , O] for Then, £(£) can be expressed as 

£(t) = (1 - R(t)W(t))Q, (15) 

where 

R(t) = l + J2(-i) n f dh f 1 dt 2 --- dt n Lse(tn)L se (t n _ 1 )---L se (t 1 ) (16) 

n=1 Jo Jo Jo 

°° ft ftl ftn-l _ 

W{t) = 1 + V i n / dh / dt 2 ■ ■ ■ / dtnQLseWQLsefo) • • • QL se (t n )Q. (17) 
~, JO Jo Jo 



and 

L se (s) ee e- <(£ " +I,e)a L a ee <(I, * +Le)i . 



These formula provide a systematic perturbative expansion. 

Finally we carry out a perturbative expansion of (1 — S(t)) -1 up to the first order of the 
interaction H se . We get the equation 

^0(t) = e lLt PtLO(0) 

- e iLt t dse- iLsS Pe- iLeS L se e i{L ° +Le)s QLO{0) 
Jo 

+ Qe iLQt iLO{0). (18) 

Notice that this is a perturbative expansion in a differential equation, so this is not a 
simple perturbation. The solution contains a certain kind of higher order contributions 
systematically in the sense of ordinary perturbation theory. 

We apply Eq.(pT8[) to a-2-n model defined by the Hamiltonian density 

h(x) = ~{$(x)$(x) + V0(x)V0(x) + m 2 s (t )<f) 2 (x)} 
+ -{U(x)U(x) + W(a;)W(aO + rr^{x)} 

+ ^0(x)vr(x)7f(x) - ^6m 2 s0 (t Q )4> 2 (x), (19) 

where 

m 2 s (t ) = m 2 s0 + 5m 2 s0 (t ), (20) 

and m S Q is the bare mass of a meson, 5m 2 (t ) is the counter term, £ is the time at which we 
renormalize the divergence, and <&(x), H(x) are the conjugate fields of <fi(x), n(x) respectively. 
Now we want to evaluate the time dependence of a meson field by regarding the ir meson 
field as environment. Thus the first, the second and the third lines on the r.h.s. of ([19|) are 
identified with H s , H e and H se , respectively. The a meson field is expanded as 

= E / 1 5 {c%{ty %S + ci(t)e-^}, (21) 
k j2Vu s (k) 



^) = ^E^ ^{^y" £ - c\(t)e-^}, (22) 



where 



u s {k) = y/fc + m? a (t ), (23) 

and V is the volume of the total system. We take the limit V — > oo at the end of the 
calculation. 

Equation of motion for the creation operator of the a meson is then given by 

jy u (t))=zu s (k)(cl(t))+zp(k) 

+ 5^§^(<4(f)> - (c_^)» 

u s {k) 

+ -^-MmA&t) - bml (t ))((cl(t)} + (c^t)}), (24) 
uJs(k) 2 

where 



p(k) = ^J= / d 3 x(7r(£,0)7r(£,0)) e e^, (25) 
pVu s $) Jv 

A(k,t) = --?- f ds [ d 3 x'd 3 x(\n(x\0)n(x',0),ii(x,s)7f(x,s)]) e 
zv Jo Jv 

xe ikx e -ikx' e -iu Js (k)s^ ^ 

7f(x,t) = e iLet n(x), (27) 

and ( ) e and ( ) mean to take expectation values by using Tr p e and Tr p, respectively. 
Here we consider the case where, at the initial time t = 0, 7r meson field is in a thermal 
free gas phase of temperature T while a meson field is in a coherent state. Therefore, the 
expectation values of 7r mesons and a mesons are taken for the thermal and the coherent 
states, respectively. The Fourier transform p(k) of tt-tt correlation function is the mean field 
felt by a mesons and is of the 1st order in g, while the 2-K-2-K correlation function A(k,t) 
is of the second order in g. The real and the imaginary parts of A(k,t) are related to a 
damping rate and a frequency shift of a meson, respectively. As seen in (|24|), equation for 
(cUt)) couples with (c_g(t)). Therefore, we need the equation for (c_g(t)): 



d 
dt 

Re A(k, t) 
u s (k) 



-A c -k(-t)) = -i^ s {k){c_^t))-ip{k) 



((#)) " <c_*(t)» 



' (W(fc,t)-^m2o(to))((4(*)) + (c_^))). (28) 



,r\ V V"! / o suv "U//W fc 

From Eg. (f26|) , we obtain 



1fi7r2|fc|yo JO 



16Ti 2 \k\ Jo Jo ^(h) s 
x {exp[i(x + uj w (ki) — u) s (k))s] + exp[—i(x + ui n (ki) + u> s (k))s] 

+ exp[i{x - uj^kt) - u> s {k))s] + exp[-i(x - u^h) + u s {k))s]}\% (29) 



where 



U = =- , (30) 



w*(k) = V^ 2 + ml, (31) 



u p = ^(\k\ + N) 2 + m 2 , (32) 



UJ, 



(|*|-|*i|) a + «& (33) 



and A is the momentum cutoff. The qualitative behavior of ImA(k, t) is shown in FIG. 1(a). 
One sees that ImA(k, t) has an ultraviolet divergence for t > 0. (KeA(k, t) which corresponds 
to the decay width of a meson is finite.) 

In order to solve the coupled equations ( f2"4|) and (|28|), one has to remove the ultraviolet 
divergence by tuning the mass counter term. We choose t = 00, that is, we require that the 
frequency shift of a meson (Aa>) vanishes at a meson momentum k = 0, T = and t = 00: 

Au = [ImA(0, 00) - ^m 2 (oo)]/w(0)| T=0 = 0. (34) 

By this procedure, the imaginary part of the renormalized 2ir-2n correlation function defined 
by lmA ren (k,t) = lmA(k,t) — |5m 2 (oo) becomes finite for t > 0, as shown in FIG. 1(b). 
The divergence of lmA(k,t) for t > is transfered to the divergence at t — after renor- 
malization. In order to see the cutoff dependence, we have carried out numerical calculation 
by taking m s = 550MeV } m n = UOMeV, g = W50MeV, k = 30MeV and T = 0, i.e., 
= 0. The result for lmA(k, t) is shown in FIG. 2. At large t, lmA(k, t) approaches rapidly 
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to A-independent limit while an appreciable A-dependence remains near t = as it diverges 
at t = as A — > oo. 

One may still worry about the presence of divergence at t — 0. However, we can show 
that it is harmless in calculating physical quantities. The reason is that the important 
quantity in solving Q2~3] ) is not A(k,t) itself but the integral Jj dsA(k, s) which is finite. In 
fact, one can show that lmA(k,t) behaves like lncrf for small but non-zero t: 

rt f A ki -> 

I ds I dk\ — [cos{u p + u^{ki) — u s (k))s 

Jo Jo suj n (ki) 

ft f A —2 

— cos(u m + uUki) — u s (k))s] ~ ds dk\ — sin(fcs)sin(2/cis) 

Jo Jo s 

. k 



/ ds- (cos2As - 1] 
Jo s 



, r sin(2At) cos (2 At) , n . 

~ H — \ — " ; . so -ln2At-7}, 35 

1 2At (2At) 2 Ul y ' 

where 7 is the Euler constant. In going from the first line to the second line on the r.h.s., 
we use the condition that k « 1/t and from the second line to the third line, At >> 1. In 
the limit of A — ► 00, we have 

the last line of (35) ~ -Mn 2At 

= -k(\n a' 1 A + In at), (36) 

where a is a constant which has a mass dimension . The first term is an ultraviolet loga- 
rithmic divergence which can be renormalized by the mass counter term. The second term 
is a new divergence which cannot be eliminated by our procedure. Fortunately, this new 
divergence is logarithmic. So the integral f dsA(k, s) is finite, and one gets a convergent 
result for (ct(t)) in the limit A — > 00. 

Furthermore it is remarkable that there is no new divergence when T 7^ if we renor- 
malize it at T = like in usual equilibrium thermal field theory. This can be seen in 
FIG. 3. 

Combining the above results, we can evaluate Re(ct(t)) for various T, of which result is 
shown in FIG. 4. Starting from an arbitrarily chosen initial value, (cAt)) makes a damped 
oscillation towards a equilibrium value which is vanishing for \k\ 7^ in this order. 



To summarize, we have calculated time dependent correlation functions using Shibata- 
Hasitsume's projection operator method. We encounter a divergence which can be renor- 
malized by a mass counter term. The remaining divergence at the initial time t = is 
harmless in calculating physical quantities. It remains as a future problem to study the 
role of wave function renormalization and vertex renormalization which will be necessary in 
more complicated cases or higher order calculations. 

Anastopoulos |15| has recently studied a similar problem by using Nakajima-Zwanzig's 
projection operator method. He uses a different projection operator and a different renor- 
malization procedure. 
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FIG.l:Qualitative behaviors of lmA(k,t) before ((a)) and after ((b)) renormalization. 



FIG.2:Cutoff dependence of lmA ren (k,t). Solid, dashed and dot-dashed lines are the 
results for momentum cutoff A = 2GeV, AGeV and 8GeV, respectively, at T = OMeV. 

FIG. 3: Temperature dependence of lmA ren (k, t) for A = 8GeV. Solid, dashed and dot- 
dashed lines are the results for temperature T = OMeV, lOOMeV and 200MeV, respectively. 

FIG.4:Temperature dependence of Re(ct(t)). Solid, dashed and dot-dashed lines are the 
results for temperature T = OMeV, lOOMeV and 200MeV , respectively, and (ct(0)) = 
IhMeV-i. 
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FIG.3 
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